Taming random lasers through active spatial control of the pump 
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Active control of the pump spatial profile is proposed to exercise control over random laser emis- 
sion. We demonstrate numerically the selection of any desired lasing mode from the emission 
spectrum. An iterative optimization method is employed, first in the regime of strong scattering 
where modes are spatially localized and can be easily selected using local pumping. Remarkably, 
this method works efficiently even in the weakly scattering regime, where strong spatial overlap of 
the modes precludes spatial selectivity. A complex optimized pump profile is found, which selects 
the desired lasing mode at the expense of others, thus demonstrating the potential of pump shaping 
for robust and controllable singlemode operation of a random laser. 

PACS numbers: 42.55.Zz,42.25.Dd 



Multiple scattering of light in random media can be 
actively manipulated through spatial shaping of the inci- 
dent wavefront[l| . This technique has allowed advances 
in focusing [2|-[J7and imaging [HI, paving the road to 
actual control of light transport in strongly scattering 
media 0-0 . Introducing gain in disordered media allows 
amplification of multiply scattered light, leading to the 
observation of random lasing The broad range of 



systems where it has been studied [ll| and the fundamen- 
tal questions it has raised 12|, [l3| has captivated the com- 
munity for this last decade. Prospective of random las- 
ing is however strongly hindered by the absence of emis- 
sion control: random lasers are highly multimode with 
unpredictable lasing frequencies and polydirectional out- 
put. Manipulation of the underlying random structure 
[1J423| and recent work constraining the range of lasing 
frequencies [3, IHI have resulted in significant progress 
toward possible control. However, the ability to choose a 
specific frequency in generic random lasing systems has 
not yet been achieved. The spatial profile of the pump 
is an interesting degree of freedom readily available in 
random laser (e.g., 26|, [27|). In a regime of very strong 
scattering where the modes of the random system are 
spatially localized [28|, local pumping allows selection 
of spatially non-overlapping modes liol , [soj . In weaker 
scattering media however (e.g., 31n23|), several hurdles 



appear toward achieving fine control. Selecting modes is 
complicated by a narrow distribution of lasing thresholds 
[3^ . and spatial mode overlap. Increased pumping 
required in these lossy systems begins to alter the ran- 
dom laser itself. Moreover, modifying the shape of the 
pump introduces changes to both the spatial and spectral 
properties of lasing modes [H, 1361-38] . Such difficulties 
are typically absent in more conventional lasers, which 
have employed pump shaping, both electrically 39l-42| 
and optically [43i], to select favorable lasing modes. The 
question is, can shaping of the incident pump field achieve 
taming of random lasers? 

In this letter, we exercise control over the distribution 



of lasing thresholds via the pump geometry to choose the 
random laser emission frequency. The fluctuation of las- 
ing thresholds from mode to mode was recently found 
to increase when pumping a smaller spatial region (i^ . 
Here we demonstrate numerically that the threshold of 
any single lasing mode for the structures investigated can 
be significantly separated from its neighbors if the spatial 
profile of the pump is correctly chosen. An iterative ap- 
proach is proposed, inspired by spatial shaping methods 
employed for coherent light control jl| in complex me- 
dia. The optimization algorithm is based on a simple 
minimization criterion and can be easily implemented 
in experiments. We first show that the algorithm sim- 
ply converges to the expected locahzed pump profile in 
the localization regime. In the weakly scattering regime, 
we show that mode selection is also possible despite the 
strong mode overlap and allows for effective monomode 
laser action. Strikingly, we find that the evolution of the 
pump profile during the optimization process alters the 
random laser in an advantageous way, so as to achieve 
the desired control. 

A one dimensional random laser is represented by a 
stack of 161 dielectric layers (optical index ni) separated 
by air gaps (no = 1.0) (see Fig.[T]). Randomness is intro- 
duced in the thickness of each layer (io,i = (<^o,i)(l + ^C)? 
where (<io,i) is the average thickness, < 77 < 1 the de- 
gree of randomness, and ^ ^ [— 1; 1] is a random number. 
The position along the system is x G [0;L], where L is 
the total length. In the following, {di) = 100 nm for the 
dielectric layers and (do) = 200 nm for air, giving a total 
average length (L) = 48.1 /im. The degree of randomness 
is set to 7^ = 0.9. 

We choose the dielectric material as host to the gain 
medium described by a frequency dependent susceptibil- 
ity 
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FIG. 1. (Color online) Principle of active control of a ID ran- 
dom laser. The black slabs represent the dielectric material 
which is also the gain medium. Optical pumping (green ar- 
rows) is transverse to the structure and amplitude modulated. 
This spatial modulation may be provided in a real experiment 
by a spatial light modulator (blue). Lasing (red arrows) oc- 
curs along the structure and will depend on the pump profile. 



density of excited atoms when the system is uniformly 
pumped, ka is the atomic transition and Aka is the spec- 
tral linewidth of the atomic resonance. As a result, the 
refractive index of the dielectric becomes complex and 
frequency dependent, ni ni{k) = ^Jn\^ Xgi^)- 
the following, the transition frequency is ka = 10.25 
/im~^ and the spectral linewidth Aka = 0.25 /im~^. The 
gain is assumed to be linear in the sense that it does 
not depend on the electrical field intensity, an assump- 
tion that is only valid at or below threshold. We use 
the transfer matrix method (sH, |46| to find the lasing fre- 
quency, the threshold and the spatial distribution of the 
lasing modes. 

When partial pumping is employed, the density of ex- 
cited atoms due to the pumping process is modulated 
according to the function /e called the pump profile as 
sketched in Fig. [TJ It can be written fE{x) x Npump^ 
where x is the position of the layer. The function fE{x) 
fulfills the constraint < fE{x) < 1 to mimic, for in- 
stance, the amplitude modulation of the pump beam by 
a spatial light modulator. This pump profile changes 
the gain provided by each dielectric layer, giving possi- 
ble control over the lasing modes of the random laser. 
Here, we aim at selecting a particular lasing mode by 
optimizing the pump profile. Experimentally, a lasing 
mode i will be selectively excited if its threshold, Ni^ is 
sufficiently low and significantly lower than that of all 
other modes. Hence, we introduce the rejection rate, 
RRi = ^'^^o^ii^j) ^ ^Jiicji compares the threshold of 
mode i with that of the mode with the lowest thresh- 
old, minj^i{Nj). Selection of mode i is achieved when 
RRi > 1 provided its threshold, Ni remains reason- 
ably low. We therefore minimize 1/RRi + aNi , with 
a properly chosen to balance each term, and apply it- 
eratively the following algorithm. A new pump profile 
fE{x) results in new lasing modes and new values of RRi 
and Ni. We apply the projected gradient method to 
fsix) 1/RRi + aNi. At each iteration, the gradi- 
ent of 1/RRi + aNi is computed. Then, fE{x) is tuned 



within [0, 1] in the direction where the projected gradient 
of 1 /RRi + aNi is minimal. Convergence is assumed if 
its relative variation is less than 10~^. 

Depending on the value of the index contrast An = 
ni — no, the random laser is studied either in strongly 
scattering regime, where light is confined well within the 
random system, or in the weakly scattering regime, where 
modes are extended. For an index contrast An = 0.60, 
we find a localization length ^ 1.7 fim <C L over the 
spectral range 10.0 fim~^ < k < 10.5 fim~^. The system 
is therefore in the localized regime and lasing modes are 
spatially confined within the system. As a rule, mode 
selection is rather easy to achieve in this case jH, [soj 
using local pumping, from the initial knowledge of the 
location of the quasimodes of the passive system. The 
localized case serves here as a test case for our iterative 
algorithm to check if modes can be selected without any 
prior knowledge of their spatial location. 

We first compute the threshold and optical frequency 
of the lasing modes for a uniform gain profile, fE{x) = 1. 
They are positioned in the frequency-threshold plane as 
crosses in Fig. [2^. Four lasing modes with reasonably 
low thresholds and partial spatial overlap are chosen 
for demonstration. Their spatial profiles are shown in 
Fig. [SJd, together with the profile of a high threshold las- 
ing mode dit k = 10.4 /im~^, associated with a lossy mode 
strongly coupled to the left end of the sample. 

We then consider Mode2 {k = 10.22 /im~^). Its re- 
jection rate for uniform pumping {fE{x) = 1) is RR2 = 
0.45 < 1, meaning it would not lase first at threshold. We 
now apply the iterative process to select this lasing mode 
{a = 10). Its rejection rate increases rapidly as shown in 
Figure [3^^. It is larger than unity after 10 iterations and 
converges to RR2 = 5.4 after 57 iterations. The relative 
increase of its threshold is less than 50%. In contrast, all 
other modes see their threshold increased by at least an 
order of magnitude. This is illustrated in Fig. [2] which 
shows the evolution of the spectrum in the complex plane. 
Mode2 has been efficiently selected. It will be the first to 
lase above threshold and the singlemode regime will be 
robust, even at relatively high pumping rate since RR2 is 
large. The optimization algorithm has been successfully 
applied to Modei^[i^3^4]. Their final rejection rates, RRi 
and threshold, Ni are given in (Hlj . 

The optimized pump profile, fE{x)^ obtained for 
Mode2 is shown in Fig. [2jd) (see also movie in ^|). It 
is similar to the lasing mode profile, g2{x) (Fig. Efc)). 
The degree of similarity is measured by the spatial cor- 
relation Cf2 = {{/e- {fE)){g2 - (^2))), where /e and g2 
have been normalized by their variance, and is close here 
to unity, C/2 = 0.82. The solution reached by the algo- 
rithm is therefore consistent with the predicted efficiency 
of a local pumping in the localized regime, even in the 
presence of moderate overlap [H, [30| . It is also interest- 
ing to note that the change of pump profile barely affects 
the frequency (Fig. [2^) and spatial profile (Fig. [2t) of the 
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FIG. 2. (Color online) Localized random laser case, (a) Fre- 
quency vs. threshold of the lasing modes at successive itera- 
tions, when optimization routine is applied to select Mode2. 
Crosses represent initial positions for uniform pumping, (b) 
Spatial profile of five localized lasing modes. The mode profile 
on the left edge corresponds to the leaky mode with the high- 
est threshold at A; = 10.4 iim~^ . (c) Spatial profile of Mode2 
after optimization, (d) Resulting optimized pump spatial pro- 
file. 



lasing modes, as expected in the localized regime [l3| . 

We have investigated in more details the working op- 
eration of the algorithm by looking at the evolution of 
the correlations Cj^^ [1^2,3,4] of the pump profile with 
Mode^^[i 2,3,4]. As the optimization routine is applied, 
Cf2 consistently increases, while correlations Cfi^^i^^,,^] 
of the other modes are progressively minimized, as shown 
in Fig [3fc, resulting in increased thresholds for these 
modes. The crosses in Fig. [3h indicate the mode with 
lowest threshold, mmj^2(^j), entering the calculation 
of RR2 at a given iteration. After working alternatively 
on Modes and Mode4, the algorithm works exclusively on 
the rejection of Mode4. This mode has the largest overlap 
with Mode2; a fine tuning of the pump profile is therefore 
required to increase its threshold without increasing too 
much Mode2's threshold. 

We turn now to the more difficult case of a weakly scat- 
tering random laser. For an index contrast of An = 0.05, 
the localization length is ^ ~ 200 /im ^ L over the fre- 
quency range 10.0 /im~^ < k < 10.5 /im~^. The system 
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FIG. 3. (Color online) (a,c) localized case, (b,d) weakly scat- 
tering case. (a,b) Rejection Rate RR2 and threshold A^2 for se- 
lected mode. (c,d) Correlations functions C/^^ [1,3,4]. Crosses 
indicates the mode with lowest threshold entering the calcu- 
lation of Ri. 



is therefore in the weakly scattering regime and lasing 
modes are spatially extended over the whole system, as 
illustrated in Fig. H^. 

Figures [3)3, d and [H present the results for the selection 
of Mode2 (k = 10.23 /im"^ and RR^^ = 0.86) (a = 0.3). 
The algorithm converges after 320 iterations with a final 
rejection rate RR2 = 2.47. Although more modest than 
in the localized case, this increase is significant enough to 
envision monomode operation of the random laser at that 
selected frequency, even when gain saturation is included 
(47| . since higher threshold modes may be suppressed. 
Similar results are obtained for all other modes tested 
(see [m ) . Figure shows the relative impact of the it- 
erative process on Mode2 and on the other modes. As 
the thresholds increase, lasing modes experience spectral 
shifts as well as spatial deformation (Fig. St), in contrast 
to the localized case. This is a direct consequence of the 
strong level of pumping required in this case and to the 
non-uniformity of the gain profile. Figure IHi shows the 
optimized pump profile (see also movie in [51|). Such a 
profile is unpredictable. A small correlation C/2 = 0.28 
is found only when comparing the stationary component 
of the mode (Fig IIJ:) [38[ and the pump profile. All the 
value of our adaptative approach is revealed here in this 
case where modes overlap is significant, threshold dis- 
tribution is fairly narrow and nonlinear effects strongly 
impact the modes. The subtlety of the optimization pro- 
cess is exemplified in Fig [3ji where the ceaseless switch 
between modes (crosses) force reduced correlation and 
even anticorrelation between the rejected modes and the 
pump profile. 

We have intentionally chosen a small index contrast 
to test the extreme case of a weakly scattering system 
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FIG. 4. (Color online) Weakly scattering random laser case, 
(a) Frequency vs. threshold of the lasing modes at succes- 
sive iterations, when optimization routine is applied to select 
Mode2. Crosses represent initial positions for uniform pump- 
ing, (b) Spatial profile of Modei (envelope) and Mode2 (full), 
(c) Stationary component of Mode2 before (envelop) and af- 
ter optimization (full), (d) Resulting optimized pump spatial 
profile. 



(^ ^ L). Intermediate regimes of scattering with var- 
ious index contrasts have been tested and mode selec- 
tion by active spatial control worked successfully as well. 
Another intentional choice is the optimization criterion 
because it is weakly constrained. Only two modes are 
compared at each iteration and as a result, rapid conver- 
gence is achieved. Stability studies of the solutions have 
been performed, which show this optimization method 
provides most often with local maxima. Better rejection 
rates may be obtained, for example, by using a global 
optimization based on the Ant algorithm [48| . We tested 
it on a a small subset of variables, which yielded promis- 
ing results. A full global optimization (over all dielectric 
layers), however, is numerically out of reach at the mo- 
ment. Though the systems studied here with a narrow 
ran ge o f emission frequencies have experimental analogs 
23, HI, typical weakly scattering random lasers yield a 
huge number of closely packed lasing modes. In such 
cases, an efficient global optimization may be necessary, 
and therefore, is under development. Finally, it is worth- 
while pointing out that the optimization problem posed 



here is different from Q since the iterative process is 
highly nonlinear : the medium itself is modified at each 
iteration by the new computed pump profile. Surpris- 
ingly, the alterations of the random laser induced by the 
pump profile are used to the advantage of the optimiza- 
tion routine, and do not hinder the goal of mode selection. 

In summary, we have demonstrated numerically that 
control of the lasing emission frequency is possible in a 
random laser. Active control of the pump profile is pro- 
posed to select any lasing mode in the emission spec- 
trum and to significantly increase the threshold of oth- 
ers. The method has been successfully tested in the 
regime of strong scattering where modes are spatially lo- 
calized and can be easily selected using local pumping. In 
the weakly scattering regime where strong spatial modal 
overlap precludes any straightforward spatial selectivity, 
the algorithm remarkably converges to a complex opti- 
mized pump profile which selects the desired lasing mode 
at the expense of the others. The proposed algorithm is 
straightforward to implement in practice. It also pro- 
vides with an optimized solution which limits the pump 
flux on the sample and, for instance, reduce optical dam- 
age. These results open further the way to active control 
of the lasing properties of random lasers by shaping of 
the pump profile. On the basis of our results, we be- 
lieve output directionality may also be achieved. For ex- 
ample, local pumping has been shown to yield unidirec- 
tional emission in random lasers through the selection of 
a new lasin g m ode (ioj generated by the presence of gain 
boundaries [50|. Additionally, relaxing the constraint of 
mode selection may allow a combination of lasing modes 
to produce the desired directionality. Moreover, because 
partial pumping modifies the spatial distributions of las- 
ing modes, the output is not bound by the parameters of 
the passive random system. For instance, the frequency 
shift observed in the weakly diffusive case could be ma- 
nipulated too to tune at will the emission frequency of 
the laser. We believe this approach will foster the inter- 
est for random lasers. But it is not restricted to this field 
and can be extended e.g. to the domain of high-power 
broad-area semiconductor lasers where it can be an alter- 
native to the issue of filamentation in order to optimize 
their brightness [Hj. 
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